Prom Finite to Infinite Range Order via Anneafing: 
The Causal Architecture of Deformation Fauhing in 
Annealed Close-Packed Crystals 
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We analyze solid-state phase transformations that occur in zinc-sulfide crystals during annealing 
using a random deformation-faulting mechanism with a very simple interaction between adjacent 
close-packed double layers. We show that, through annealing, infinite-range structures emerge from 
initially short-range crystal order. That is, widely separated layers carry structurally significant 
information and so layer stacking cannot be completely described by any finite-range Markov process. 
We compare our results to two experimental diffraction spectra, finding excellent agreement. 

PACS numbers: 61.72.Dd, 61.10.Nz, 61.43.-j, 81.30.Hd 
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There has been considerable interest in understand- 
ing planar disorder in crystals for some iiiae^'^'^'^'^'^ 
There are several reasons for this emphasis. Since pla- 
nar defects correspond to a shift of an entire layer of 
atoms that nonetheless preserves the crystallinity within 
the layer, the disorder is confined along the stacking di- 
rection and results in a structure that may be treated as 
one-dimensional, making theoretical analysis tractable.^ 
Additionally, different stacking configurations can affect 
the physical properties of the material. For example, 
the band gap in SiC changes as the stacking structure is 
changed® and an anomalous photovoltaic effect has been 
observed in ZnS crystals that have disordered stacking 
sequences. 

While many crystals exhibit planar defects, they are 
especially common in a class of materials known as poly- 
types!^ Polytypism is the phenomenon where a crystal is 
built up from the stacking of identical two-dimensional 
layers, called modular layers (MLs)iSiiiii2ii^ Typically 
there is a small, finite set of possible ways a ML may be 
stacked upon another and, since these different stacking 
orientations often preserve the atomic coordination num- 
ber for nearest and next-nearest neighbors, the energy 
difference between different stackings can be quite small. 
It is perhaps not surprising then that polytypic materi- 
als have many different crystalline structures. There are, 
for instance, 185 known crystalline stacking structures in 
ZnS, some with unit cells extending over 100 MLsi2ii4 
Other highly polytypic materials include SiC, Cdl2, and 
Agl with about 150, 200, and 50 known crystalline struc- 
tures each, respectivelyiS 

Understanding the variety and origin of spatial orga- 
nization in crystalline ZnS polytypes on length scales 
clearly in excess of the calculated inter-ML interactions — 
- 1 ML for ZnS^s and - 3 MLs for SiC^fe— has been a 
puzzle for some time and numerous theories have been 
proposedji^iiLiSiiS Recently, we extended this by demon- 
strating that disordered ZnS polytypes also possess a 
long-range spatial organization in excess of the calculated 
inter-ML interactioniiiii^ 



In this paper we simulate transformations between or- 
dered and disordered polytypes, in particular, proposing 
a simple model to analyze the solid-state transformation 
of annealed ZnS crystals from a 2H to a 3C structure and 
finding a novel structural description — in the form of an 
e-niachinc^"'^^ — of the resulting disordered twinned 3C 
crystal. Our model is decidedly simpler than previous 
ones, designed with the goal of determining the minimal 
complication necessary to produce the experimentally ob- 
served long-range spatial order. For example, from the 
e-machine we show that infinite-range spatial memories 
arise even though the interaction length between MLs is 
restricted to nearest neighbors. 

Disordered crystals are often formed by stressing a per- 
fect crystal thermally, mechanically, or through irradia- 
tion so that stacking faults are introduced. Various mod- 
els have been proposed to explain these transformations. 
Typically, the faulting process is assumed to proceed 
slowly, so that a ML is randomly chosen and a fault is 
introduced if the local stacking configuration meets some 
criteria, usually derived from a (theoretically or empir- 
ically determined) Hamiltonian that describes the ener- 
getics of the inter-ML interactions. Further, a particular 
faulting mechanism is assumed. 

Many authors have previously performed simulation 
studies on transformations in polytypes to understand 
the structure of faulted crystals. For example, by as- 
suming interactions between MLs of up to a distance 
of three, Kabra and Pandey^^ were able to deduce that 
layer-displacement faulting was the primary defect found 
in the 211 6H^^ transformation in SiC. They also dis- 
covered that this faulting mechanism can result in long- 
range correlations between MLs without short-range or- 
der. Engel24 applied a similar model to treat the 2H 
— » 3C transformation in ZnS via deformation faulting, 
with the assumption of next-nearest neighbor interac- 
tions. These authors were able find stacking configura- 
tions that gave x-ray diffractograms qualitatively similar 
to experimental ones. Shrestha et alm^ treated the 2H 
—>■ 3C martensitic transformation in hypothetical close- 
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packed structures by introducing the infinitely strong re- 
pulsion model. They assumed that the presence of a 
fault inhibited the introduction of another fault on ad- 
jacent MLs. That is, that there was some coordination 
between faults, called nonrandom faulting. '^^•'^'^ Gosk^-^* 
developed a model to describe disordered 2H crystals that 
used a probability function for faulting on the next ML 
as one scans the crystal that depended on the distance 
from the last observed fault. By tuning various model 
parameters one could introduce nonrandom faulting into 
the stacking sequence. 

ANNNI modelsi^ have also been used to explain poly- 
typism since they are known to give long-range spatial or- 
ganization with only short-range interactions (up to next- 
nearest neighbor). However, ANNNI models require fine 
tuning the coupling parameters in an interaction Hamil- 
tonian and, thus, may be of limited applicability to poly- 
types. 

ZnS has two stable phases: the hexagonal close-packed 
(or 2H)i24i^ for temperatures above 1024 C and the cubic 
closed-packed (or 3C) for temperatures below 1024 Ci^^ If 
2H ZnS crystals are cooled sufficiently fast to room tem- 
perature, they will retain this 2H structure even though 
it is not the stable phase; presumably due to the difficulty 
of shifting MLs to form the cubic (3C) structures. These 
crystals are often further annealed for an hour at temper- 
atures 500 — 800 C, thus perhaps providing the necessary 
activation energy to induce a solid-state transformation 
from the 2H to a twinned 3C structure. This transforma- 
tion has been studied before, and many authors have con- 
cluded that deformation faulting is the main mechanism 
responsible fiSiSS at least for the case of weak faulting. 

Engel and Needsi^ performed a first-principles pseu- 
dopotential energy calculation at T = to determine the 
coupling constants between MLs up to separation three. 
They found an expression for the energy in the form 

2 

m— 1 2—1 

where the Si G {0, 1} are inter-ML spinal and Jm 
the inter-ML coupling constants at separation m. They 
found Ji > and J2 < 0, but much smaller in magnitude 
than Ji . All other couplings between MLs were found to 
be negligible. 

Our approach is to assume the following very simple 
model describes the 2H — > 3C transformation in ZnS. 

(i) The energetics of ML stacking in ZnS is describable 
by a one- dimensional Ising chain with a nearest-neighbor 
interaction Hamiltonian Ti. between MLs. Further, we 
assume that entropic effects are small at all temperatures 
of interest, so that the free energy F = H. 

(ii) The effective coupling Ji(T) between adjacent MLs 
is temperature dependent, being positive at temperatures 
below the transition, Tc ~ 1000 C, and negative above. 
While little is known about the temperature dependence 
of the effective couplings Ji{T)^ this assumption gives 



TABLE I: The deformation-faulting update rule for ZnS crys- 
tals suggested by the Hamiltonian Ti with Ji{T) > and 
■hiT) ~ 0. Only those spin flips that lower the energy of 
the spin configuration are allowed. This is equivalent to ap- 
plying elementary cellular automaton (EGA) rule 232 to the 
spin neighborhood r)\ = (s*_i, s*, s*_|_i). s\^^ gives the value 
assigned to the center spin when it is selected for updating. 
The rule is applied asynchronously to randomly selected sites. 



nt 


111 110 101 100 on 010 001 000 




1110 10 



the correct stable phases for ZnS (211 and 3C) above and 
below the the transition temperature. 

(iii) There is only one faulting mechanism — deformation 
faulting — and it is driven by Glauber dynamics. Defor- 
mation faulting^^ occurs in close-packed crystals when 
there is a slip by a non-Bravais lattice vector between 
adjacent MLs. In terms of the spin sequence, this cor- 
responds to flipping a single spin^^ and so the faulting 
reduces to basic Glauber dynamics.''*' 

(iv) The transformation is slow, sluggish, and random. 
By slow, we mean that the time for a slip between MLs 
to occur is much shorter than the interval between slips. 
By sluggish, we mean that a slip will only occur if it is 
energetically favorable. And by random, we mean that 
each ML is equally likely to be selected for possible fault- 
ing, without any coordination between MLs. Since a de- 
formation fault occurs when one portion of the crystal 
slips relative to the other, it is reasonable to assume that 
this happens infrequently and only when the slip results 
in a reduction of the energy. While the mechanism of 
faulting assumed here is expressly random, we will show 
nonetheless that this leads to a crystal structure that has 
a nonrandom fault distribution. 

Our model is most similar to that introduced by En- 
geli24 He assumed two kinds of faults, fast and slow, de- 
pending on the configuration of the five-spin neighbor- 
hood centered about the candidate spin to be flipped. 
Restricting the range of interaction to nearest-neighbors, 
as we do, only has the effect of suppressing the slow fault 
mechanism. 

We simulate the solid-state transformation in ZnS by 
assuming an initial configuration of a pure 2H crystal 
with N = 1048577 MLs, which is quickly brought into a 
temperature range (500 — 800 C) so that the 3C structure 
becomes the stable phase, but there is sufficient thermal 
energy for MLs to slip. We randomly choose a spin in the 
configuration and apply the update rule listed in Table 
We visit the spins in a round-robin fashion, so that each 
is updated once and only once, until all have been. It 
is easy to see that the update rule permits a spin to flip 
only once, if at all, and so we lose no generality in this 
model by visiting each spin only once. 

Let us define the faulting parameter f as the ratio of 
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FIG. 1: Configurations during the 2H — > 3C transforma- 
tion under asyncfironous celiular automaton 232 on a 100-site 
lattice with periodic boundary conditions. The vertical axis 
shows the degree / of faulting and the horizontal axis is the 
spatial configuration, with black squares indicating Si = 1 
inter-ML spins and white squares Si = inter-ML spins. Ini- 
tially, the crystal is in the 2H configuration as indicated by 
the alternating black and white squares at the top of the dia- 
gram (/ = 0) . The faulting proceeds in a round-robin fashion 
until all sites are visited once. At the bottom of the dia- 
gram, (/ = 1), the crystal is fully faulted and there are only 
odd-spin domains (see text) with each domain having at least 
three spins. 

number Uapp of times the update rule has been applied 
to a crystal of size N: f — Uapp/N. The faulting pa- 
rameter provides a useful index of the amount of faulting 
the crystal has been subjected to or, equivalently, of the 
stage of the transformation. That is, one can determine 
the crystal structure at various stages along the transfor- 
mation by selecting different values for / G [0,1]. This 
algorithm is formally equivalent to applying elementary 
cellular automaton rule 232 asynchronously to a fraction 
/ of randomly chosen spins (but no spin visited more 
than once). We then call this model asynchronous cellu- 
lar automaton 232 (ACA 232). Note that the update rule 
is spin-inversion symmetric, in accord with the Hamilto- 
nian. Models of this basic flavor have been used to study 
systems other than polytypism, such as cluster growth 
and phase separation in one dimension ^^i^ diffusion- 
reaction problems in one dimension^Si^ and the voter 
model. '^'^ The transformation under ACA 232 of spatial 
configurations as a function of / is shown in Fig. ^ start- 
ing from the period-two (2H) configuration. 

We then extract a structural description of a spin 
configuration in the form of an e-machine, as a func- 
tion of /, using the causal-state splitting reconstruction 
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FIG. 2: The reconstructed e-machine that describes the archi- 
tecture of the stacking process of a faulted 2H crystal across 
the range of possible fault parameters / £ [0, 1]. The nodes 
represent causal states and the directed arcs, transitions be- 
tween them. The edge labels s\p indicate that a transition 
occurs from one causal state to another on symbol s with 
probability p. 



(CSSR)^ algorithm. Specifically, we considered crystals 
with various amounts of faulting / e [0,1] in increments 
of A/ — 0.10 and performed e-machine reconstruction 
for each /. We limit the number of causal states in the 
reconstructed e-machine by requiring that the addition 
of a new causal state reduce the entropy rate by at least 
1%. 

The e-machine for a faulted 2H crystal appears in 
Fig. [3 It has 10 causal states and is a function of 3 
transition probabilities — pi , p2 , and ps . These transition 
probabilities are themselves a function of the faulting pa- 
rameter / and this dependence is shown in Fig. |31 For 
small values of the faulting parameter, pi w /. This 
implies that for small /, the e-machine largely cycles 
between the causal states A and B. This produces the 
spin sequence ...101010..., which we recognize as the 
2H stacking sequence. So, for small faulting, the origi- 
nal 2H crystal remains largely intact; as one expects. At 
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FIG. 3: The causal-state transition probabilities pi,P2, and 
P3 (see Fig. I^J as a function of the faulting parameter / 
for a crystal {N = 1048577 MLs) initially in a periodic 2H 
stacking sequence. The Pi axe estimated from reconstructed 
e- machines. 



the other extreme, when / = 1 (a fully faulted crystal) 
we have pi = 1, and there is no longer a transition be- 
tween the causal states A and B. Since there are no other 
causal-state cycles^^ that generate the . . . 101010 . . . spin 
sequence, we see that the original 2H structure is elimi- 
nated. 

Since the e-machine has spin-inversion symmetry — i.e., 
it is invariant under 1 <^ exchange, we only need to ex- 
amine half of it to get an intuitive understanding of the 
structures it captures. Let us consider only the causal 
states A, B, D, F, I, and J. As stated before, the causal 
state cycle [AB] generates the 2H structure. At B, how- 
ever, it is possible to emit a and make a transition to 
D. This implies that the spin history upon entering D is 
100 and, conversely, having seen the spin sequence 100 
uniquely sets the e-machine in D. At this point, we see 
that D must make a transition to F on spin 0. That is, 
the sequence 1001 does not occur. At F, the stacking 
process may either enter A on a 1 giving the cumula- 
tive spin sequence 10001 or may proceed to I on a 0. 
At I, the process must see another and transition to 
J, giving the cumulative spin sequence of 100000. This 
implies that the spin sequence 100001 likewise does not 
occur. In fact, we see that the causal states D and I 
force generating a as the next spin and, thus, prevent 
sequences where two or four Os are sandwiched between 
two Is. Similar reasoning shows that any even number of 
Os (Is) cannot be trapped between two Is (Os). In other 
words, the sets of sequences {10^"1, n = 1, 2, 3, . . .} and 
{01^"0, n = 1,2,3,.. .}, and any configurations that con- 
tain them, are forbidden. Thus, starting from a short- 
range crystalline structure (211 or . . . 010101 . . .), anneal- 
ing has produced complex spatial structures of infinite 
range. 

We can understand how annealing accomplished this 
increase in structural complexity in the following way. 
Define a spin domain as a sequence of k consecutive like 



spins sandwiched between two unlike spins, where fc is a 
positive integer. In the original 2H structure, there are 
no even- A; spin domains. The update rule has the effect of 
joining two spin domains of the same spin by eliminating 
the k — I spin domain separating them. Therefore, the 
resulting domain is of size knew — /c; + /c,. + 1, where fc; 
(kr) is the k value of the spin domain to the left (right) 
of the flipped spin. The spin that is flipped also con- 
tributes a site to the new spin domain. Thus, if ki and 
kr are initially odd, knew is also. Since an even spin do- 
main cannot appear, as one scans the crystal one must 
remember the number of consecutive like spins in order to 
determine the admissibility of the next spin value. In gen- 
eral, this can require remembering an indefinite number 
of previous spins. In this way, the process has an infinite 
memory, so we say that the memory length r is infinite. 
Put another way, spins arbitrarily far apart may contain 
information about each other not contained by the in- 
tervening spins. Symbolic dynamics calls such a process 
strictly soficm^ The most significant consequence of being 
sofic is that no finite-order Markov process can describe 
the stacking sequence. Notably, this suggests that previ- 
ous efforts to describe annealed, disordered crystals that 
assume some sort of Markov niodel2i2i^iS*ii*iSiii22i2& are 
likely inadequate. 

As we noted earlier, the structure for a completely 
transformed crystal occurs for / = 1. Then the short- 
est causal-state cycle is [FACEBD] and generates the spin 
configuration . . . 111000 . . .. There can be forays into the 
causal states I and J or G and H, but these only serve to 
increase the spin domains by an even number. Thus, the 
final crystal is a disordered, twinned 3C crystal with only 
odd spin domains. Notably, some time ago, Mardij*^ 
performed a statistical analysis of the known crystalline 
polytypes in ZnS and found a significant bias toward odd 
number Zhdanov elements in the configurations. These 
correspond to regions of odd spin domains, in agreement 
with odd spin domains predicted by ACA 232 for disor- 
dered ZnS crystals. 

It is useful to inquire into the origin of the infinite 
memory length. Clearly it does not derive directly from 
the interaction Hamiltonian, which has an interaction 
range of one. In fact, it comes from the repeated in- 
teraction of the stochastic (asynchronous) dynamics and 
the structural constraints imposed by the allowed (de- 
formation) faulting mechanism. When a spin is flipped, 
information is lost. That is, when scanning the spin con- 
figuration, the process must remember the last 01 or 10 
pair. The local update rule acts to eliminate pairs of this 
form and so requires the process to remember an increas- 
ingly longer "history" of spins, as annealing progresses. 

We now address the issue of nonrandom faulting. By 
examining the update rule in Tabled we see that when a 
spin flips in the ACA 232 model, the two adjacent spins 
are necessarily prevented from also flipping. This there- 
fore precludes the possibility of two adjacent faults. This 
implies that faults are not random, since there must be at 
least a two ML separation between them. This is easily 
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FIG. 4: The diffraction patterns along the 10. Z row for par- 
tially annealed 2H crystals. All spectra, here and elsewhere 
(except for those from the fault model), are normalized to 
unity. 



seen by examining a sequence of spins and the conditional 
probabilities for faulting as implied by the e-machine. 

Take a segment of an unfaulted 2H crystal as 10101010. 
Now introduce a fault (i.e. spin flip) in the spin sequence 
to get 10001010, where the underline indicates the flipped 
spin. Now neither spin adjacent to the flipped spin can 
itself flip. However, the spins at a distance two or greater 
can. Let us consider additional spin flips at a separation 
of two, three, and four. We then have the sequences 
10000010, 10001110, and 10001000 as possible stacking 
sequences, representing deformation stacking faults at 
separations of two, three, and four. As before, the his- 
tory 100 requires the e-machine to be in causal state D. 
The next spin must be 0, advancing the e-machine to 
causal state F. Thus the causal state architecture of the 
e-machine has prevent two adjacent faults, as required, 
and the spin sequence is now 1000. If the stacking se- 
quence is to resume unfaulted, we expect the next spin 
to be 1. From the e-machine, however, the probability of 
0, implying a deformation fault separated from the first 
by two MLs, is p2- That is, the conditional probabil- 
ity of a second fault occurring two MLs after the first is 
P2- Similar reasoning shows the e-machine specifies that 
the conditional probability of a deformation fault occur- 
ring three and four MLs after the first is pi. We find, 
therefore, that the conditional probabilities of observing 
deformation faults at a separation of one, two, and three 
MLs are all different, even though AC A 232 manifestly 
assumes only interactions between nearest neighbors. For 
small faulting, however, Fig. |31 shows that pi ~ P2, so 
that this distinction is minor for faults at a separation of 
two and three. 

The above exercise demonstrates explicitly that the e- 
machine can represent a nonrandom distribution of stack- 
ing faults and this distribution can arise from an inter- 
action range of one combined with a simple mechanism 
for annealing. ACA 232 then provides some theoreti- 
cal justification for models that assume structures with 



a nonrandom fault distribution without appealing to any 
long-range interactions or coordination between faults. 

To see how well ACA 232 describes phase transitions 
in experimentally observed in ZnS crystals, we calculated 
the diffraction spectrum along the 10. Z row for the faulted 
crystal for various values of the fault parameter as shown 
in Fig. ^ (We employ the same procedure, definitions, 
and assumptions here as we used previouslyi^) Let I 
be a dimensionless variable that indexes the magnitude 
of the perpendicular component of the diffracted wave. 
For small values of the faulting, the 211 Bragg peaks at 
1 = 1/2 and I = 1 are widened considerably. As / ap- 
proaches 0.50, the enhancement in diffracted intensity at 
these positions has vanished and there is instead signifi- 
cantly increased scattering near the Z w 1/3 and Z « 2/3 
positions, which is normally associated with twinned 3C 
structure. There are also small rises in the diffracted in- 
tensity at the / « 1/6 and Z « 5/6 positions, which is 
often considered a sign of 6H structure. As / approaches 
1.0, the peaks at I « 1/3 and I « 2/3 sharpen and the 
diffracted intensity at the Z « 1/6 and Z « 5/6 posi- 
tions weaken. For / = 1, the crystal has found a local 
minimum of the Hamiltonian, and no further transfor- 
mation occurs. The crystal is now a disordered, twinned 
3C structure. 

Notice also that it had been assumed previously that 
the small rises in the intensity at I ~ 1/6 and Z w 5/6 
signaled that parts of the crystal were faulting into the 
6H structure, but we see that no special assumption of 
this kind is needed. Rather, the spin sequence 000111 
(i.e. the 611 structure) occurs quite naturally as the 2H- 
crystal deformation faults, giving the enhanced diffracted 
intensity at these positions. 

We now compare the predictions from our simple 
model with two experimental ZnS diffraction spectra we 
have previously treatedii*i^: the first from an unannealed 
sample; the second, from an annealed sample. Figure [S] 
shows the spectrum from an as-grown, disordered 2H 
crystal. Our previous analysis showed that, while there 
were displacement faults, growth faults, and some 3C 
structure, the dominant faulting was due to deformation. 
It is reasonable, then, to ask if ACA 232 can reproduce 
this spectrum. For / = 0.06 ACA 232 does reproduce the 
Bragg-like peaks at I « 1/2 and / « 1 rather well. The 
small rise in diffuse scattering at Z « 2/3 is not recovered, 
nor should one expect it to be. This feature comes from 
the small amount of 3C structure present and, as we saw 
above, ACA 232 is not able to reproduce this. That it 
does not lends further strength to our earlier interpre- 
tation^^ that a small amount of 3C structure is present. 
In fact, we find that for small faulting, ACA 232 gives 
nearly the same results as the fault model, as can be seen 
in Fig.El 

Figure |^ compares an experimental ZnS diffraction 
spectrum from a 211 crystal annealed at 500 C for 1 h with 
an / = 85% deformation-faulted 211 simulated diffraction 
spectra. The agreement is excellent. The deformation- 
faulted crystal reproduces the two disordered 3C peaks 
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FIG. 5: A comparison of an experimental diffractogrami— *ii 
along the 10. Z row (triangles) with the fault model (solid 
line>ii and ACA 232 (dotted line) at / = 0.06. Except for 
some minor differences, ACA 232 gives the same description 
as the fault model. We find the profile 7?.-factos— between 
experiment and ACA 232 to be 29%. We see that ACA 232 
describes the Bragg peaks at I « 0.5 and I « 1 well, but misses 
the small rise in intensity at I « 0.67. This likely results from 
3C structure, so it is not surprising that it is absent the ACA 
232 spectrum at low /. 



in the experimental spectrum a,t I w —2/3 and / « —1/3, 
as well as the diffuse scattering between these diffrac- 
tion maxima. Even the slight enhancement in diffraction 
at / ~ —1/6 and I w 1/6 is captured (see inset). Our 
deformation-faulting model gives a stacking configura- 
tion that is spin symmetric and so it is unable to account 
for the relative difference in intensity between the two 
peaks at I « —2/3 and I « —1/3. There is a question here 
concerning the quahty of the data,^^ and so the difference 
in intensities may be an experimental artifact. Note that 
the diffracted intensity between the two Bragg-like peaks 
at I « —2/3 and I « —1/3 is especially well represented 
by ACA 232. This contrasts with the fault model which 
simply cannot recover the diffuse scattering. 

There are however, features seen in diffractograms of 
partially annealed ZnS crystals that are not observed 
in the simulated diffractograms of Fig. 0| Namely, the 
diffractograms of some partially transformed ZnS crys- 
tals, as in Fig. 9 of Ref^^- , show enhancement simulta- 
neously at values of I associated with the 2H structure, 
namely Z « 1/2 and Z « 1, as well as at / values associ- 
ated with twinned 3C structure (/ ~ 1/3 and I ~ 2/3). 
In our simulation studies, the strong 2H reflections at 
I 1/2 and I « 1 disappear before the 3C reflections at 
Z w 1/3 and Z w 2/3 start becoming prominent. ACA 232 
therefore has difficulty even qualitatively explaining such 
spectra and one is led to entertain models more compli- 
cated than ACA 232. 

We introduced a very simple model, called ACA 232 
due to its equivalence with elementary cellular automa- 
ton rule 232, to study the 211 3C phase transition 
in polytypic ZnS crystals. We assumed only nearest- 



FIG. 6: A comparison of an experimental diffractogram-*— 
along the 10./ row (triangles) with the fault model (solid 
line)2*ii and ACA 232 (dashed line) at / = 0.85. With a 
profile 7^-factor of 18% (compare 13% for eMSR^ and 33% 
for the fault model), ACA 232 captures a surprisingly large 
amount of the stacking structure. Notice particularly that the 
inset shows that ACA 232 nicely reproduces the small rise in 
diffracted intensity at I « —0.16. 



neighbor interactions between MLs and thus a local, 
though asynchronously applied, faulting mechanism. We 
simulated a solid-state transition during annealing of the 
crystal and found a description of the stacking structure 
for all amounts of faulting in the form of an e-machine. 
The resulting e-machine implied that, despite interac- 
tions only between nearest neighbors, the structure of 
the crystal had a nonrandom distribution of faults. The 
e-machine also shows an infinite memory length, which 
means that arbitrarily separated spins share nonredun- 
dant information. The corresponding symbolic dynam- 
ics is strictly sofic and this implies that no finite-order 
Markov process can completely specify the stacking in 
these faulted crystals. Intuitively, this arises because 
only odd-length spin domains are allowed, if the crystal 
starts from a perfect 2H specimen. The infinite mem- 
ory length process discovered here contrasts with the 
r — 3 memory length we found by analyzing disordered 
ZnS crystals using 

eMSR.i3 This suggests that while the 
r = 3 e-machine captures much of the stacking struc- 
ture present, it is nonetheless only an approximation to 
a more structurally complex stacking process. The re- 
lation between the e-machine found in these simulation 
studies and that determined by eMSR from experimen- 
tal spectra is currently a subject of investigation. Since 
even simple, annealed crystals probably have an infinite 
memory length, this further suggests that the Markov 
approximations used in the eMSR algorithm should be 
circumvented and a more general spectral reconstruction 
technique is warranted. This, too, is a current subject of 
research. Finally, we compared diffraction spectra calcu- 
lated from ACA 232 to previously analyzed ZnS spectra 
and found good quantitative agreement. 

While not included here, our studies also considered 
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different initial stacking sequences. We found that of- 
ten the memory length was infinite after annealing and 
that the structure of the final faulted crystal was highly 
dependent on the initial stacking structure. Although 
phase transitions in actual ZnS crystals undoubtedly are 
more complicated, our simple model captures much of 
the noncritical organization they undergo. Also, consid- 
ering its simplicity, and the fact that similar models have 
been used to study other systems, the results presented 
here should have applicability far beyond the realm of 
polytypism. 
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